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The development of advanced quantum technologies and the quest for a deeper understanding of 
many-particle quantum mechanics requires control over the quantum state of interacting particles to 
a high degree of fidelity. However, the quickly increasing density of the spectrum, together with the 
appearance of crossings in time-dependent processes, makes any effort to control the system hard and 
resource intensive. Here we show that in trapped systems regimes can exist, in which isolated energy 
bands appear that allow to easily generalize known single-particle techniques. We demonstrate this 
for the well-known spatial adiabatic passage effect, which can control the center-of-mass state of 
atoms with high fidelity. 
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I. INTRODUCTION 

Understanding the effects of interactions between 
many particles at the quantum level is an important 
task for increasing the access to, and control over, ever 
larger parts of the Hilbert space mm. However, this is a 
difficult problem, as interacting many-particle states are 
usually too complex to allow for exact analytical treat¬ 
ment and often require numerical resources that are out 
of reach for current computers. One way to approach this 
problem is to study small systems first and use the de¬ 
veloped understanding for scaling up. This allows theo¬ 
retical treatment and experimental verification, as recent 
progress in cold atomic gases has led to control over the 
trapping of small numbers of particles and the ability to 
measure them with high precision [3]. Using Feshbach 
resonances one can then change the interaction strength 
between the particles and engineer well-defined quantum 
states. The importance of this lies in the ability to ex¬ 
plore new aspects of fundamental quantum mechanics as 
well as to design new applications in quantum technolo¬ 
gies. 

However, no generic and straightforward strategies ex¬ 
ist to develop new engineering tools that take advantage 
of the full Hilbert space. One controlled way to make 
progress is to try to generalize known techniques for 
single-particle states to either weakly-correlated many- 
body states or to small strongly-correlated systems !3]. 
In this work we focus on the well-known spatial adia¬ 
batic passage (SAP) protocol 0-0, which is a technique 
that allows to transfer a localized wave function between 
different positions in space by adiabatically following a 
specific energy eigenfunction. SAP for ultracold atoms 
has up to now only been studied for a single atom 003- 
lllj . weakly-interacting systems H2HE], or fermionized 
bosons m- Here we discuss a system of two interacting 
bosons in a triple well potential and, using exact solu¬ 
tions and Hubbard models, we show that the interactions 
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can destroy, but also revive, the possibility for SAP due 
to an interaction-induced energy-band separation. This 
separation allows to generalize single-particle processes 
to multiple particles and requires to consider repulsively- 
bound pair co-tunnelling processes Q2HSI- 

In the following we will first very briefly review the 
technique of spatial adiabatic passage (Section 0 and 
the general solution for two interacting particles in a har¬ 
monic trap (Section III). We then numerically solve the 
problem of SAP for interacting particles (Section |lV[ ) and 
discuss the limiting cases for weak and strong interactions 
using the Hubbard model (Sections |V| and [vl| . The im¬ 
portance of level crossings and non-adiabatic evolution is 
discussed in Section IVIII before we conclude. 


II. SINGLE-PARTICLE SAP 

The SAP protocol for a single atom involves three de¬ 
generate localized trapping states | j) with j = l, m and 
r (for left, middle and right), which are centered at po¬ 
sitions d,L < d,M < d-R. The system is described by the 
Hamiltonian 

Ho = LM\l)(m\ + n M R\m)(r\ + h.c. , (1) 

where the time-dependent nearest-neighbor couplings 
Cljji are controled by the distance between the traps [5j. 
Note that throughout the paper we use dimensionless 
units where h, atomic masses, and trapping frequencies 
are equal to 1. One of the eigenstates of H 0 is the so- 
called dark state [25] . 


\D) = cos9\l) — sin0|r), (2) 

with tan$ = UIlm/^mr and SAP describes the trans¬ 
port of a particle from |Z) to |r) following | D) by chang¬ 
ing 9 from 0 ( UL M r > H L m ) to 7r/2 (f 1 M r < H L m )• 
This is achieved with the trap movement shown in 
Fig- 0a), which maintains an energy gap between the 
d ark state and the two other eigenstates of the order of 
V^lm + n MR ( see Fig- 0 b)) and changes the dark state 
from 1 1) at t = 0 to |r) at final time t = T (see Fig. 0c)). 






2 



0 0.2 0.4 0.6 0.8 1 

t/T 


FIG. 1. (Color online) (a) Positions of the three harmonic 
well minima for the SAP protocol as used in our simulations 
(dashed red: dn, dotted green: cIm, dot-dashed blue: dij). 
The initial (and final) distance between wells is d max = 9, the 
minimum distance is d m i n = 3 and the time delay between 
the two approaches is T/10. (b) Energy eigenvalues of the 
single-particle Hamiltonian IT]), with the one corresponding 
to | D) displayed in blue (gray), (c) Coefficients of | D) in the 
{|jr)} basis (dashed red: |Z), dotted green: | m), dot-dashed 
blue: |r)). 


To avoid excitations and ensure that SAP succeeds, the 
whole process needs to be carried out adiabatically. 


III. EXACT MODEL FOR TWO PARTICLES 


For simplicity, we will use a one-dimensional model 
of two interacting bosons in a triple-well potential, even 
though our results can easily be extended to higher di¬ 
mensions. The Hamiltonian is then given by 

H = +V(x k ,t)\ +gS(x!-x 2 ), (3) 

k =1 ' k ' 


where Xk is the position of the fc-th atom. The trapping 
potential V is modeled by a piecewise harmonic triple 
well 02 with minima located at d l, dM, and d r (see 
Fig. ma)). The last term describes the contact inter¬ 
action (of strength g ) between the atoms [2B]. 

From here on, we refer to the ground state of two inter¬ 
acting bosons in each harmonic well as | j) (j = L, M, R ), 
whose energy E g is related to g by [ZB) 


2V2T(1 - E g /2) 
F((l - E g )/ 2) 


(4) 


where r(:r) is the gamma function. Throughout the pa¬ 
per we will discuss the whole range of E g : from the 
non-interacting case (g = 0, E g = 1) to the Tonks- 
Girardeau (TG) limit (g —> oo, E g = 2). Our objective 
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FIG. 2. (Color online) Final population in state \R) as a 
function of E g after the two-particle SAP protocol is carried 
out over a total time T = 4000 (blue/dark gray) or T = 12000 
(orange/light gray). Dotted vertical lines indicate energies for 
which the spectrum is shown in Figs. [3] and [4] 


is to find an eigenstate of H, analogous to the single¬ 
particle dark state (tending to | L) and |i?) at initial and 
final times) which allows for two-particle SAP transport 
between these states. 


IV. SAP FOR INTERACTING PARTICLES 


With the system initially in state | L), we numerically 
integrate the time-dependent Schrodinger equation with 
the Hamiltonian ([3]) while changing the positions of the 
traps as shown in Fig. [lja). We then calculate the SAP 
fidelity F (population of |i?) at final time t = T) for E g 
between 1 and 2. The results for two different values of 
T are shown in Fig. [2] where one can identify several 
regions of interest. 

The process succeeds for the non-interacting case 
{E g = 1), where the particles are independent and each 
undergoes single-particle SAP. Full transfer also occurs 
in the TG limit {E g = 2), which can be easily under¬ 
stood since the bosonic atoms can be treated as indepen¬ 
dent fermions occupying the two lowest single-particle 
energy levels, and each of them can then be treated with 
a Hamiltonian similar to Eq. 0 El- The fidelity drops 
sharply for energies near these extreme values, regions A 
and C in Fig. [2j where the population transfer is only 
partial and depends on T. However, from E g ~ 1.12 to 
E g ~ 1.45 (region B) a plateau appears where F > 0.998. 

In the rest of this work, we will analyze these three re¬ 
gions in detail by diagonalizing the exact Hamiltonian of 
the system. We will also introduce two Hubbard Hamil¬ 
tonians, for weak and strong interactions, whose formal¬ 
ism will allow us to study the structure of the eigen¬ 
states and give additional insight into the importance of 
co-tunneling processes. 
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V. WEAK INTERACTIONS 


Restricting the analysis to two-particle states in the 
lowest Bloch band, the system can be modeled for weak 
interactions using a finite-size Bose-Hubbard Hamilto¬ 
nian [27] 


h b = £ 

j—L,M,R 


u 

—n,- 
2 J 


( n j - 1 ) + egrij 


(5) 


+ £Ilm (b\b M + b' M b L ) + Q.mr (b\fbR + 


q( c °) 

T “LM 


/jt 2 /) 2 _i_ /,t 2 *,2 \ I q( co ) (hf 2 h 2 h* 2 h 
°L °M + °M°L) + “Mi? [°M°R + °R 0 


where foj and bj are the creation and annihilation op¬ 
erators for a boson in the ground state of well j and 
rij = b^bj is the number operator. The onsite interaction 
is described by XJ and the ground state energy is eo = 1/2. 
Single-particle tunneling rates (fijy) and two-particle co¬ 
tunneling rates ) between wells j and j' are numer¬ 
ically calculated via the Gram-Schmidt orthonormaliza¬ 
tion procedure [9j using the single-particle [5. and two- 
particle [26] wave functions, respectively. 

When the two traps are close to each other, Q n < and 

are of the same order of magnitude, but for weak 
interactions single-particle tunneling dominates due to 
its larger amplitude. For stronger interactions, however, 
single-particle tunneling becomes off-resonant, unlike co¬ 
tunneling which dominates because it remains resonant. 
While co-tunneling is not usually considered, it is crucial 
to understand the behavior observed in the previous sec¬ 
tion. A study of Hb (with real-time simulations, eigen¬ 
values and eigenstates) with and without co-tunneling 
terms can be found in the appendix. 

In order to understand regions A and B in Fig. [2] we 
diagonalize numerically both H and Hg at any point dur¬ 
ing the SAP evolution. The spectra of H for E g = 1.05 
and 1.25 are shown in Fig. [3](a,b), and consist of three 
distinct bands. The lowest band contains three states 
(where the two atoms sit in different wells), and for large 
trap separations (i.e., initial and final times) has energy 
1. Not discussed here, this band contains an eigenstate 
which allows the adiabatic transfer of an atomic hole 
between the outermost wells for strong enough interac¬ 
tions nu. The middle band (also with three states) is the 
one of interest to us because it only involves states with 
the two atoms in the same trap, \j), with energies around 
1 + U = Eg. The higher band, with energies around 2, 
consists of six states which correspond to states where 
the atoms sit in the ground state and first excited states 
of different traps. For the same parameters, the spec¬ 
trum of Hb is shown in Fig. [3](c,d) . These spectra show 
two bands, corresponding (with good agreement) to the 
two lower bands of the exact Hamiltonian we have just 
discussed. Higher bands do not appear because only the 
lowest Bloch band has been considered in Hg. 

The band around E g has a similar structure to the 
single-particle SAP spectrum (cf. Fig. 0b)) and it 



FIG. 3. (Color online) (a,b) Lowest 12 eigenvalues of H for 
the two-particle SAP scheme with the trap moving sequence 
of Fig. 0a) for (a) E g = 1.05 and (b) E g = 1.25. In (b) 
the energy of the dark state (asymptotically | L) and |i?)) and 
the state with which it couples the most are marked in blue 
(dark gray) and orange (light gray), respectively. The inset 
shows a zoom-in of the marked crossing between these two 
states (marked with a circle) which is analyzed in Sec. |VII| 
(c,d) Eigenvalues of Hb for the same parameters as (a,b), 
with the energy of the dark state drawn in blue (gray). (e,f) 
Coefficients in the Fock basis |{nj}) of the dark states in 
(c,d) (dashed red corresponds to \L), dotted green to | M), 
dot-dashed blue to | R), and the solid lines to states where the 
two atoms are in different traps). 


contains an eigenstate, drawn in blue (dark gray) in 
Fig. §c ,d), which allows for SAP transport because of 
its dark-state-like structure, shown in Fig. [3je,f ). It is 
important to remark that this state is only present if co¬ 
tunneling is considered in Hg, see the appendix. The 
SAP transport fails for weak interactions, Fig. [3])a,c,e), 
because when tunneling becomes relevant the two bands 
overlap and dynamically following the dark state is hard 
due to the presence of level crossings. We can then see 
that the appearance of the plateau in Fig. [2] for stronger 
interactions, i.e., U > y/2fl ~ 0.15, is due to the two 
bands remaining separated during the whole process, see 
Fig- 0b ,d). This allows for the SAP process to be suc¬ 
cessful again because the band with the dark state re¬ 
sembles the single-particle SAP spectrum, cf. Fig. Qb). 
The dark state has a similar shape to the single-particle 
one but for the states | j) coupled through co-tunneling 
(compare Figs. [3](f ) and 0c)). 
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VI. STRONG INTERACTIONS 

One expects then that for stronger interactions the pro¬ 
cess will keep working, until E g approaches the band with 
energies around 2, which will cause the fidelity to drop 
again in region C of Fig. [2] Because of the fermionic 
behavior of the TG gas, this regime is best modeled by 
restricting the system to two-particle Fock states with 
one atom in each of the two lowest Bloch bands and us¬ 
ing a finite-size Fermi-Hubbard Hamiltonian m 


Hf= £ 

j=L,M,R 


Urijoriji + ^2 


i=0,l 


( 6 ) 


^LM a Li a Mi + ^MR a \li a Ri + h-C. 


i=0,l 


+ Q\ J M®L0 a \,l a M0 a Ml + h.C. 

, rv( co ) ft i i 

E ^MR a MO a Ml a RO a Rl + h.C. 


Here a;„- and a,-,- are the fermionic creation and annilii- 

jx jx 

lation operators for a particle at energy level i of well j, 
iljj, are the tunneling frequencies at level i, riji = 
and ti = 1/2 + i. The onsite interaction is given by U, 
which is now negative since two bosons interacting repul¬ 
sively with finite strength have less energy than two non¬ 
interacting fermions. In other words, the ground state 
energy is now e 0 + ei + U = 2 — \U\ = E g . 

We show in Fig.|4^a,b) the spectrum of H for E g = 1.6 
and 1.85, and one can see the same three-band struc¬ 
ture as in Fig. [3] For the same values of E g , the spec¬ 
trum of Hr, Fig. [4jc,d), consists of two bands: around 
2 — \U\ = E g and 2, in good agreement with the exact 
results. Once again, one state in the band around E g (in 
blue/gray) has the form of a dark state, see Fig. 

As E g increases, the two bands start to overlap and cross¬ 
ings appear, leading to the dark state involving states of 
the upper band, see Fig. |4](b,d) , which perturb the adia¬ 
batic transfer as seen in region C. It is worth noting that 
Hr contains a dark state even without the co-tunneling 
term. 


VII. (A)DIABATIC REQUIREMENTS 

Thus far we have seen that our system’s Hamiltonian 
contains a band with energies around E g which consists of 
a three-level system of states | j), coupled by repulsively- 
bound pair tunneling. Thanks to the atomic interaction, 
this band can be isolated from the rest of two-particle 
states, and allows the transport to succeed. From the 
size of the two bands of Hr, one would expect that pro¬ 
cess should not be affected by level crossings until around 
E g ~ 1.7. However, the plateau in Fig. [2] only extends 
until E g ~ 1.5. In order to understand this we will look 
back on the spectrum of the exact Hamiltonian spectra in 
Figs. §b) and[4ja). For both spectra we have highlighted 



FIG. 4. (Color online) (a,b) Lowest 12 eigenvalues of H for 
the two-particle SAP scheme with the trap moving sequence 
of Fig. [I|a) for (a) E g = 1.6 and (b) E g = 1.85. In (a) the 
energy of the dark state (asymptotically | L) and |R)) and 
the state with which it couples the most are marked in blue 
(dark gray) and orange (light gray), respectively. The inset 
shows a zoom-in of the marked crossing between these two 
states (marked with a circle) which is analyzed in Sec. |VII| 
(c,d) Eigenvalues of Hr for the same parameters as (a,b), 
with the energy of the dark state drawn in blue (gray). (e,f) 
Coefficients of the dark state in (c,d) the Fock basis \{n g i}) 
(color coding is the same as in Figs. Hfe,f)). 


in blue (dark gray) the energy of the dark state, i.e., the 
state which at initial and final times is | L) and | R), re¬ 
spectively, and therefore the state we are initially follow¬ 
ing. In contrast to the Hubbard Hamiltonians, this dark 
state crosses another eigenstate (shown in orange/light 
gray) twice, creating a finite probability for the system 
to leave the dark state. For 1.1 < E g < 1.8 these are the 
only relevant crossings affecting the dark state, and we 
examine them carefully in the following. However, as the 
crossings are actually avoided crossings, it is clear that 
the speed at which they are passed will determine the 
adiabaticity and the success of the transfer [28} [29]. 

The probability for a system following an eigenvector 
|i(f)) to be excited into another one \j(t)) between times 
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FIG. 5. (Color online) Probability of transition at the 
crossing between the eigenstates shown in the insets of Figs. [3] 
and [4] for different total times T and energies E g . Dashed 
vertical lines indicate the total times used in Fig. [2] 


to and tf can be approximated by [30] 


Pi-tj — 


Ito 0‘(*)l 

±\i{t))e r *o {E ^- E ^ dT dt 


Ito 

2 


2 


(7) 


where E/.(t) is the energy of state | k(t)). The denomina¬ 
tor was added from the original expression for normaliza¬ 
tion purposes and represents the probability for the state 
to be excited for an infinitely fast process. Full trans¬ 
fer through the SAP protocol can be achieved if both 
crossings are passed either adiabatically {pi^j = 0, fol¬ 
lowing always the dark state) or completely diabatically 
= 1, following the orange (light gray) state between 
the crossings). For intermediate values of Pi~>j , the sys¬ 
tem’s actual state will be distributed between different 
eigenstates and the transfer will not be complete. 

In Fig. [5] we show the calculated transition probabili¬ 
ties for the gap as a function of T and E g . One can see 
that for the timescales of T used in Fig. [2] the trans¬ 
fer is completely diabatic for E g = 1.25. For E g = 1.6, 
however, Pi^j starts to falls to 0.99 (0.96) for T = 4000 
(12000), and even lower for higher E g . Therefore, the 
plateau ends because of the increased size of the gap in 
the avoided crossing, which no longer allows for diabatic 
passage. Moreover, it can be seen that the length and po¬ 
sition of the plateau is T-dependent, since the transition 
from a diabatic to an adiabatic process happens at very 
different time scales, which depend on E g . For longer 
time scales, e.g. of the order of 10 7 , two plateaus would 
appear in region B: one from E g ~ 1.1-1.4 in which the 
transfer in this crossing would be completely diabatic, 
and one at E g ~ 1.5-1.7 where the transfer would be 
adiabatic. 
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FIG. 6. (Color online) Populations of the eigenstates from 
integrating the Schrodinger equation with the Hamiltonian 
Hb for Eg = 1.25 (a) without and (b) with the co-tunneling 
terms. The initial state is \ip(t = 0)) = | L) and the tunnel¬ 
ing rates are calculated numerically with the time-dependent 
trap positions in Fig. [lja). Color coding is the same as in 
Figs.[3}e,f). 



FIG. 7. (Color online) Spectra of (a) the exact Hamiltonian 
(lowest six eigenstates), and Hb without (b) and with (c) the 
co-tunneling terms for the SAP process as a function of time 
for Eg = 1.25 (U = 0.25). The energy of the dark state in (b) 
is shown in blue (gray). 


continuous region for intermediate interactions exist over 
which high fidelities can be obtained. This is due to 
the fact that for intermediate values of E g a decoupled 
energy band appears, which possess a dark state facili¬ 
tated by two-particle co-tunneling. However, when this 
band overlaps with other energy bands, the appearance of 
level crossings prevents the robust use of the dark state. 
This behavior is generic to any multi-well setting and not 
specific to SAP. It is worth noting that the above effect 
is limited to systems where no phonon modes exist and 
therefore does, for example, not apply to SAP in quan¬ 
tum dot systems. 


ACKNOWLEDGMENTS 


This work was supported by the Okinawa Institute of 
Science and Technology Graduate University. We thank 
Irina Reshodko for helpful discussions. 


VIII. CONCLUSION 

In this work we have studied the spatial adiabatic pas¬ 
sage protocol for a system of interacting bosons over the 
entire range of repulsive interactions. We have found 
that, in addition to the trivial cases for non-interacting 
and infinitely strongly interacting particles, a large and 


APPENDIX: JUSTIFICATION OF THE 
CO-TUNNELING TERMS 

In this appendix we will examine the validity of the 
Bose-Hubbard Hamiltonian in greater detail and con¬ 
firm the importance of the co-tunneling terms. For this 
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FIG. 8. (Color online) Populations of the eigenstates of the 
Hb without the co-tunneling terms in the Fock basis for the 
SAP process as a function of time. States are ordered with 
increasing energy as (a) to (f), and use the same color coding 
as Figs-IJe.f). 
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FIG. 9. (Color online) Same as Fig. [8] but including the co¬ 
tunneling terms in Hb- 


we have performed numerical simulations of the time- 


dependent Schrodinger equation using Hb for E g = 1.25 
(U = 0.25) and, because this falls inside region B 
of Fig. [2] the process is expected to work. The re¬ 
sults, shown in Fig. [6j compare the cases where the co¬ 
tunneling terms in Eq. (|5j) are either neglected or consid¬ 
ered and one can immediately notice that the simulation 
that does not consider the co-tunneling terms gives in¬ 
complete transfer to state |i?). However, when one con¬ 
siders the co-tunneling terms, the transfer is complete, 
which establishes that these terms are crucial in order to 
explain the two-particle dynamics. 

To shed more light on this we have also computed the 
spectrum of Hb , without and with co-tunneling, and that 
of the exact Hamiltonian during SAP. The results are 
shown in Fig. [TJ and one can see that the inclusion of the 
co-tunneling terms leads to Hb approximating the exact 
spectrum more closely. Furthermore, we have computed 
the populations of the eigenstates in the Fock basis with¬ 
out (Fig. [8]) and with (Fig. [9]) co-tunneling terms. States 
in the lower band, depicted in (a-c), are mostly com¬ 
posed by states where the atoms are in separate traps. 
Because these state do not couple through co-tunneling, 
both their structure and energies coincide in Figs. [8] and 
[9] It is noteworthy that the state in (b) corresponds to 
another kind of dark state that transfers an atomic hole 
from the right well to the left well [15] . The states in 
the higher band, shown in (d-f), clearly differ for the two 
Hamiltonians. When the co-tunneling terms are absent, 
no dark state that allows transitions from | L) to |7?) ex¬ 
ists in Hb, see Fig. [8] However, it is clearly present when 
taking into account co-tunneling, see Fig. [9](e) (the en¬ 
ergy of the dark state is shown in blue (gray) in Fig.[7^b)). 

Therefore we have clearly established that a dark state 
that allows transfer between |L) and |I?) is only present 
in Hb if the co-tunneling terms are considered. 
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